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Abstract — A new method for removing point radio sources and other non-Gaussian noise is proposed as a means 
of improving the accuracy of estimates of the angular power spectrum of the cosmic microwave background (CMB). 
The main idea of the method is to reconstruct fluctuations of the CMB in places contaminated by such emission, while 
traditional methods simply exclude these regions from consideration, leading to the appearance of "holes"in the resulting 
maps. The fundamental possibility of reconstructing the CMB signal in such holes follows from the analytical properties 
of a function with a finite spatial spectrum (the Silk damping frequency). A two-dimensional median filter is used to 
localize the point radio sources.Results of simulations of the method for maps of modest size are presented. The efficiency 
of applying the method to reconstruct the CMB from data with limited resolution and contaminated by appreciable pixel 
noise is investigated. The fundamental possibility of applying the method to reconstruct the CMB distribution in the 
region of the Galaxy is also demonstrated. @2005 Pleiades Publishing,Ine. 


1. INTRODUCTION 

In connection with numerous experiments on mea¬ 
suring the anisotropy of the cosmic microwave back¬ 
ground (CMB) with the aim of deriving high-accuracy 
estimates of the main cosmological parameters of the 
Universe [1], it is currently of interest to find new ap¬ 
proaches to data reduction that could increase the ac¬ 
curacy with which various background contaminating 
signals of Galactic and extragalactic origin can be re¬ 
moved from radio observations of the CMB.This prob¬ 
lem has been widely discussed in the literature [2-4]. 

We will restrict our consideration to the removal of 
pointlike radio sources and other non-Gaussian noise, 
which most strongly distorts the angular power spec¬ 
trum of the CMB at high spatial harmonics (multi¬ 
poles). This problem has been considered by a number 
of authors [5-9].The currently available methods differ 
primarily in the means used to localize the contami¬ 
nating signal, with the aim of subsequently removing 
the corresponding sections from the measured maps. 

The following three strategies for the removal of 
pointlike sources are all natural: (1) reconstructing the 
point sources,then subtracting them from the measured 
CMB maps; (2) excluding contaminated pixels from 
the maps; and (3) reconstructing the true values of the 
CMB in contaminated pixels. 

Experience shows that the first strategy leads to 
larger errors in the angular power spectrum of the CMB 


than the other two methods, since it is not possible to 
completely accurately reconstruct the brightness dis¬ 
tribution of the point sources when the spatial spectra 
of the individual distinguished signals overlap [10]. For 
this reason, maps obtained with this method are char¬ 
acterized by the presence of appreciable non-Gaussian 
residuals. The second approach, which requires only 
the determination of the locations of contaminated pix¬ 
els and their elimination before estimating the CMB 
power spectrum, is substantially more effective. Pre¬ 
cisely this method is currently the most widely used. 

However, this second method for the removal of 
contaminating point sources is likewise not free of 
drawbacks. First, the total effective area of the map 
that is used to estimate the angular power spectrum is 
reduced, leading to appreciable errors for more distant 
sections of the CMB. Second, the simple elimination of 
contaminated pixels does not fully restore the Gaus¬ 
sian statistics of the CMB, since it leaves ’’holes”in the 
maps that are used. In addition, the sharp boundaries 
of the cutout sections lead to so-called Gibbs phenom¬ 
ena in the behavior of the spectrum, which are espe¬ 
cially strongly manifest at high spatial harmonics,and 
can appreciably complicate studies of the secondary 
anisotropy of the CMB. 

Note as well that increasing the resolving power of a 
system used to measure the CMB anisotropy will lead 
to the detection of a larger number of weak sources 
[11], so that an increasingly large number of pixels in 




the CMB map are contaminated by these signals. A 
special problem is also presented by the strong con¬ 
tamination due to Galactic emission. Simple removal 
of the zone of Galactic emission, which comprises an 
appreciable fraction of the overall area of the celestial 
sky, unavoidably leads to appreciable errors when es¬ 
timating the angular power spectrum compared to the 
use of completely uncontaminated data. 

Obviously, if the third strategy is realizable in prac¬ 
tice, which depends primarily on the noise character¬ 
istics of the apparatus, it will be free of the backs of 
the first two methods.Therefore,we propose and devel¬ 
op here a new method for the removal of high-multipole 
non-Gaussian noise, including both the localization of 
such noise and reconstruction of the true values of 
the CMB anisotropy in the contaminated regions. For 
the specific conditions considered in this paper, this 
could lead to an increase in the accuracy with which 
the angular spectrum of the CMB can be construct¬ 
ed,compared to the simple elimination of contaminated 
map pixels. 

Thus,our aim is to present and develop a new 
method for eliminating the contamination of a CMB 
signal by point sources and other sources of high- 
multipole noise. We estimate the internal accuracy of 
the method and its stability to the input noise, and also 
investigate its application to data obtained with finite 
resolution for various levels of the instrumental pixel 
noise, thereby evaluating possibilities for applying the 
method to real systems. 

The following sections of the paper present a statis¬ 
tical description of the CMB anisotropy; describe our 
model for the observed map, the method used to local¬ 
ize the point sources, and the method used to recon¬ 
struct the CMB values in noise-contaminated regions; 
and present the results of numerical simulations car¬ 
ried out for maps with modest angular dimensions but 
having the same statistical properties as the CMB dis¬ 
tribution over the entire celestial sphere. 

2. MODEL FOR AN OBSERVED MAP 
OF THE CMB 

The distribution of the CMB temperature over the 
celestial sphere can be presented by the following ex¬ 
pansion in spherical harmonics 17 ^( 6 >, 0 ) [12]: 

AT 
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where T and AT are the mean temperature and tem¬ 
perature fluctuations of the CMB and the a^rn are the 
coefficients of the expansion. The angular power spec¬ 
trum of the fluctuations is determined as the mean 
square of the coefficients airn’- 
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where i is the multipole number. 


If the fluctuations in the early Universe satisfy 
Gaussian statistics, as is expected in most cosmolog¬ 
ical theories, each coefficient aira should be statistical¬ 
ly independent. Thus, the power spectrum Cg provides 
a full statistical description of the CMB anisotropy, 
which is a fundamental characteristic of the Universe 
that can be obtained directly from observations via a 
spherical-harmonic analysis. 

For Gaussian fields, the expansion coefficients 
^ 7^ On likewise represent Gaussian fields with ran¬ 
dom phases, zero means,and the dispersions 

< aimCL^/^f >= ^ii' ^rnrn'Cl^ i ^ 0. 

Here, we consider sections of the sky that are mod¬ 
est in size. In this case, it is expedient to operate with 
the Gaussian fields in a flat, two-dimensional space, 
and to replace the spherical-harmonic analysis with 
a Fourier analysis, which appreciably simplifies the 
testing of proposed methods. The CMB fluctuations 
AT = T— < T > can then be generated via the simple 
calculation of the Fourier series [12] 
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where L is the linear size of the region considered in 
radians, (Ox^Oy) are Cartesian coordinates on the sky 
(in the spatial region), and is the number of 

the Fourier component in the region with 

spatial frequencies u and v. 

The amplitudes of the Fourier components 
D{nu^ riy) obey a Gaussian distribution with zero mean 
and the dispersion 

< \D{nu,ny)\'^ >= Ci, ^ + n2, (2) 

while the phases are uniformly distributed in the inter¬ 
val (0,27r). Here, Ci is the angular power spectrum of 
the CMB temperature when it is expanded in spherical 
harmonics. 

Relation (2) describes the circular symmetry of the 
power spectrum, i.e., its independence of the azimuthal 
number m, which we will use below as additional a pri¬ 
ori information when reconstructing the CMB fluctu¬ 
ations in noise-contaminated regions. 

Another, more important, characteristic of the pow¬ 
er spectrum is its finite spatial extent,which follows 
from the existence of the so-called Silk damping fre¬ 
quency above which the power-spectrum fluc¬ 

tuations fall off sharply and the contribution of the 
CMB to the total observed signal becomes negligibly 
small with increasing frequency. The finiteness of the 
spectrum enables us to apply the theory of analytical 



functions [14] to describe the CMB fluctuations, which 
implies the possibility of reconstructing functions over 
an entire region based on knowledge for part of the re¬ 
gion or at some set of points [15]. This principle lies 
at the basis of the proposed algorithm for reconstruct¬ 
ing the CMB fluctuations in regions contaminated by 
noise. 

When modeling the background of point sources 
(PS),we assume that they are randomly distributed 
over the sky in accordance with a Poisson law.The in¬ 
strumental pixel noise is ’’white’’Gaussian noise with 
zero mean.The measured CMB map can be represent¬ 
ed by the model 

CMBmeas = {CMB + PS) * BEAM -h N, (3) 

where N represents noise, BEAM the antenna beam, 
and * a linear convolution. 

The required data processing consists of solving an 
equation of the form of the convolution (3) for the CMB 
signal. The sequence of operations we have used during 
the reconstruction of the CMB is the following. 

(1) The noise is filtered using a Wiener filter or a 
modification of such a filter that does not distort the 
form of the angular power spectrum [16]; this yields an 
estimate of the signal {CMB -f- PS) * BEAM. 

(2) This quantity is deconvolved from the antenna 
beam using a regularized inversion filter [17,18], in or¬ 
der to derive an estimate of the signal CMB + PS. 

(3) Positions contaminated by point sources are lo¬ 
calized using a median filter (see Section 4). 

(4) The CMB signal is reconstructed at the contam¬ 
inated positions using the method proposed in the fol¬ 
lowing section. 

3. RECONSTRUCTION 
OF THE CMB FLUCTUATIONS 
AT CONTAMINATED POSITIONS 

The proposed method for reconstructing the fluctu¬ 
ations in holes in a CMB map is a modification of the 
algorithm of Fienup [19], intended for the reconstruc¬ 
tion of an image of an object with a finite carrier from 
the amplitude of its Fourier spectra (the phase prob¬ 
lem). In our modified version, we determine the limi¬ 
tation on the spectral region using information about 
the finiteness and circular symmetry (2) of the ampli¬ 
tude spectrum,rather than the amplitude spectrum it¬ 
self, and determine the limitation on the spatial region 
using the known map values. 

The algorithm is iterative, and consists of the fol¬ 
lowing sequence of operations. 

(1) An initial approximation for the map is made. 
We recommend the use of the initial map with zero 
brightness in specified locations (holes) as a first ap¬ 
proximation. 


(2) The Fourier transform of the initial approxima¬ 
tion is calculated, bringing about a transformation to 
the spatial-frequency domain. 

(3) A condition for the spatial limitation of the fluc¬ 
tuation spectra,which comes about because the values 
of the Fourier components derived in the first step with 
numbers £ > £d have been set equal to zero,is im¬ 
posed.An additional constraint on the region of spatial 
frequencies that appreciably speeds up the convergence 
of the algorithm is the circular symmetry of the power 
spectrum of the fluctuations. To satisfy this last con¬ 
dition, the form of the Fourier components is modified 
so that the power spectrum has a form that is consis¬ 
tent with relation (2). For each value of £, the squared 
amplitudes of the spectra measurements are averaged 
in radius over a length £, after which these measure¬ 
ments are replaced by values with the derived squared 
amplitudes, retaining the phases. 

(4) The inverse Fourier transform of the spectrum 
obtained in the previous step is taken, bringing about 
a transformation to the domain of the CMB map. 

(5) The constraints on the spatial region of the map 
are imposed. The brightness values outside the holes in 
the map are replaced by the known values, while the 
values inside holes are not changed. 

(6) The Fourier transform of the map obtained in 
step 5 is taken. 

(7) Return to step 3 until the image obtained in 
step 5 ceases to change in accordance with a chosen 
convergence criterion. 

Simulations show (see Section 6) that this modified 
Fienup algorithm leads rather rapidly to the correct so¬ 
lution. For this to be the case, the number of unknown 
values of the CMB fluctuations should be approximate¬ 
ly half the number of knowns, if the discretization fre¬ 
quency of the map is twice its upper spatial frequen¬ 
cy. If the discretization frequency is increased, a larger 
number of map values can be reconstructed, since the 
insufficient information in the spatial domain can be 
compensated by information in the spatial-frequency 
domain. 

Obviously, a high computational speed for the 
method can be reached via the application of Fast 
Fourier Transform algorithms. 

4. METHOD FOR LOCATING 
CONTAMINATED SECTIONS 

The proposed method for localizing the point 
sources is based on applying a two-dimensional median 
filtration,which is a row-by-row column algorithm con¬ 
sisting of one-dimensional n-point median filtrations 
[20,21]. The number n is taken to be odd. If n = 2/c + I, 
one-dimensional median filtration consists of ascribing 
to the current value for the sequence the mean of the 
series that is obtained when the (2/c-|-I)-point sequence 



is placed in increasing order, with the first k values lo¬ 
cated to the left and the last k values to the right of 
the given value. As a result of applying this operation 
to a one-dimensional series, all impulsive noise is elim¬ 
inated. In our case, such noise corresponds to point 
sources. 

To obtain the series of removed sources itself, we 
subtract the output map from the input map. Obvious¬ 
ly, since the nonlinear transform is applied to the useful 
signal together with the noise, the resulting difference 
map obtained from the subtraction contains additional 
noise, whose magnitude is appreciably lower than the 
desired impulsive noise when n is chosen correctly. Sec¬ 
tions of the resulting difference map that are contam¬ 
inated by point sources can be identified by applying 
a cutoff to the levels of this map. The lower the cutoff 
level, the more pixels are subject to distortion. In our 
case, the cutoff level was usually chosen to be fairly 
low (10% of the peak value of the CMB fluctuations), 
in order to avoid missing genuine places that are con¬ 
taminated by weak radio sources. On the other hand, 
the number of identified points must not be too large 
to ensure convergence of the subsequent reconstruction 
of the CMB-fluctuation values at these points. 

In contrast to numerous methods that have been 
proposed in the literature (wavelet analysis [7,8], clean 
and maximum-entropy deconvolution [10], optimal lin¬ 
ear filtration [9], etc.), the proposed method for localiz¬ 
ing the point sources is distinguished by its simplicity 
and the speed with which it operates, which is impor¬ 
tant in the reduction of large datasets obtained over 
the entire celestial sphere. 

Obviously, this method will operate more reliably 
the brighter the point sources and the fewer the num¬ 
ber of point sources in the map being analyzed. We 
have investigated the method for the cases of both 
relatively bright (resolved) and relatively weak (un¬ 
resolved) radio sources, whose number grows with the 
resolution of the instrument used [11]. 

5. FILTRATION 

OF THE INSTRUMENTAL NOISE 
AND DECONVOLUTION 

We used a modification of a Wiener filter that pre¬ 
serves the form of the angular power spectrum [16], 
called a power filter, for the filtration of the instru¬ 
mental white noise. A power filter can be applied in 
two ways. In the first, the dispersion of the instrumen¬ 
tal (pixel) white noise is taken to be known a priori^ 
whereas in the second, it is not. We present a prelim¬ 
inary estimate of this noise based on applying a me¬ 
dian filter to the input signal. The dispersion of the 
instrumental noise was determined by subtracting the 
output signal of the median filter from the measured 
map, CMBmeas- We present here results obtained us¬ 
ing the second approach. Reconstruction of the map 
taking into account the antenna beam was carried out 


using a regularized inversion filter containing a reg¬ 
ularizing parameter that depends on the level of the 
residuals after filtration of the noise [17,18]. 

6. SIMULATION RESULTS 

In this section, we present the results of applying 
the proposed method to model CMB measurements. 
The anisotropy of the CMB was simulated using a func¬ 
tion for the angular power spectrum that corresponds 
to a standard ACDM cosmological model for the Uni¬ 
verse with D5h^=0.02, Da=0.65, D^=0.3, h=0.65, and 
n=l. We numerically generated a map of the CMB 
7.5*^ X 7.5*^ in size in accordance with formulas (l)and 
(2); this map, therefore, possesses all the properties of 
the CMB fluctuations distributed over the entire celes¬ 
tial sphere. The upper frequency of the angular pow¬ 
er spectrum of the CMB corresponds to the multipole 
^=1536. 

6.1. Effect on the CMB Reconstruction 
in Point-Source Contaminated Locations 
in Systems with High Angular Resolution 

The goal of this experiment (experiment 1) is to il¬ 
lustrate the proposed method and to estimate the effect 
of applying the method in systems with high angular 
resolution that detect a large number of relatively weak 
radio sources.We tested the method assuming that the 
preliminary noise filtration and deconvolution were al¬ 
ready carried out (see the end of Section 2). This was 
done with the aim of estimating the internal accuracy 
of specifically the method for reconstructing the CMB 
in holes in the initial CMB map. 

We took the map to be 64x64 pixels in size, which 
corresponds to discretization at the Nyquist frequen¬ 
cy for the specified maximum multipole of the CMB 
spectrum,^^aa^ = 1536. The results of this numerical 
experiment are shown in Fig.l, which presents (a) 
the initial (model)map of the CMB fluctuations; (b) 
the map of the CMB contaminated by radio sources 
{CMB + PS] 234 weak sources with amplitudes ex¬ 
ceeding the CMB-fluctuation peak by a factor of three, 
distributed randomly over the map according to a Pois¬ 
son distribution); (c) the map obtained by processing 
the map in (b) with a row-column, three-point medi¬ 
an filter; (d) the map obtained by replacing the CMB 
values at the detected contaminated points with ze¬ 
ro brightness; (e) identified sections of contamination 
(which, in general, differ from the intrinsic distribution 
of point sources, see Section 3)obtained by applying an 
amplitude cutoff of 10% of the peak CMB fluctuation 
to the difference of maps (b) and (c); and (f) the map 
of the CMB reconstructed using the modified Fienup 
algorithm,which nearly precisely coincides with the in¬ 
put model CMB map in (a). 

Figure 2 presents the angular power spectra of the 
maps; the solid curve refers to the original and recon¬ 
structed CMB maps, the dotted curve to the contam- 



mated map (b), the short-dashed curve to the map 
obtained via median filtration (c), and the long-dashed 
curve to the map with zero brightness in contaminated 
locations (d). The dispersions of the signals represent 
a quantitative characteristic of the experiment,and are 
presented in Table 1. We can estimate the effect of 
applying our approach compared to the traditional ap¬ 
proach that does not reconstruct the CMB fluctuations 
in holes by comparing the angular power spectra in 
Fig.2 (as well as the dispersions of the corresponding 
residual maps). 

We can see that simple removal of 234 contaminat¬ 
ed measurements from the CMB map leads to distor¬ 
tion of the true angular power spectrum, especially at 
high harmonics: this distortion is from 30% to 100% 
when ^>1000. Applying the proposed method to re¬ 
move the point sources leads to virtually no such dis¬ 
tortion. Simulation of the case of a background of unre¬ 
solved point sources with amplitudes a factor of three 
smaller than the maximum CMB fluctuations likewise 
yielded very high accuracy for the reconstruction in the 
absence of pixel noise. 

To investigate the stability of the method with re¬ 
gard to the input instrumental noise,we added Gaus¬ 
sian white noise with a signal-to-noise ratio of SNR^^IO 
to the map in panel (b), taking the signal level to be 
that of the CMB fluctuations. The dash-dotted curve 
in Fig.2 depicts the angular power spectrum of the re¬ 
constructed map in this case. We can see that small 
variations in the input data gave rise to only small 
changes in the solution; i.e., the proposed algorithm 
displays a fairly high stability. Analysis of these results 
indicates that simple zeroing of the CMB values, even 
in precisely determined contaminated locations, yields 
appreciable errors in the reconstruction of the angular 
power spectrum, especially on small angular scales. If 
the contaminated sections occupy a large area (Fig.le), 
the error in the desired angular power spectrum can be 
still larger, possibly leading to substantial errors in de¬ 
rived estimates of cosmological parameters. 

Thus, in this experiment, without allowance for the 
real resolution of a system and the input instrumental 
noise, it was possible to reach very high accuracy in 
the CMB-map reconstruction. As we indicated above, 
this testifies to a high internal accuracy of the pro¬ 
posed method for reconstructing the CMB in holes in 
a map. Note that convergence of the algorithm to the 
required solution can also be achieved without plac¬ 
ing constraints on the circular symmetry of the power- 
spectrum fluctuations.However,use of this information 
as an additional constraint appreciably speeds up the 
process of convergence. Experience shows that,when 
the map is discretized with a frequency that is twice 
the bandwidth of the spectrum signal, high accuracy 
in the reconstruction is achieved only if the number of 
unknowns is no more than one third of the total num¬ 
ber of pixels in the map. As is shown above,the method 


displays a fairly high degree of stability to the noise in 
the input data; i.e., small variations in these data cor¬ 
respond to only small variations in the reconstruction. 
When the input noise level is high, a preliminary filtra¬ 
tion should be applied using a power filter (see Section 
5), as was done in the experiments discussed in the 
following section. 

6.2. Effect on the CMB Reconstruction 

in Contaminated Locations in the Model 

CMBmeas=(CMB + PS )^BEAM + N 

Let us make the problem more complex and con¬ 
sider a model for the observed CMB map that cor¬ 
responds to expression (3), which includes the effect 
of the limited resolution of the system and the input 
pixel noise.In the following set of experiments,the map 
size was 128x128 pixels, the pixel size was 3.51', and 
the full width at half maximum of the antenna beam 
10.53'. For an antenna with a diameter of 1 m, this 
resolution corresponds to an observing frequency of 98 
GHz.These system parameters are close to those for 
one of the channels of the PLANCK mission [11]. 

We carried out four experiments based on these pa¬ 
rameters (experiment group 2), with various levels of 
white pixel noise,beginning with zero noise. The dis¬ 
persions of the input,noise,reconstructed,and residual 
maps can be used as a measure of the accuracy of the 
reconstruction of the CMB anisotropy, and are giv¬ 
en in Table 2. The results are shown in Fig.3, where 
the numbers (columns)indicate the experiment number 
and the letters (rows) indicate the physical meaning 
of the map: (a) measured CMB map satisfying (3);( 
b) reconstructed CMB map with holes in contaminat¬ 
ed pixels, including the effect of the antenna beam; 
(c) residual map, equal to the difference of the given 
CMB map and the reconstructed map shown in row 
”b”; (d) reconstructed CMB map obtained using our 
method to interpolate functions in the holes; (e) resid¬ 
ual map,equal to the difference of the given CMB map 
and the reconstructed map shown in row ”d”. Obvious¬ 
ly, each column of maps presented in the figure corre¬ 
sponds to a particular value for the pixel-noise level. 

The results presented in Fig.3 were obtained using 
pixel-noise filtration based on applying a median filter 
to derive a preliminary estimate of the noise dispersion 
(see Section 5). The effect of the antenna beam in the 
reconstructed map was included using a regularized in¬ 
version filter. 

Analysis of the results presented in Fig.3 shows that 
the main advantage of the proposed method for remov¬ 
ing point sources is that the residual maps shown in 
row ”e”are free from non-Gaussian features, and resem¬ 
ble residual Gaussian noise, whose amplitude depends 
on the level of the input Gaussian noise (see the dis¬ 
persion in Table 2). As we can see in map row ”c”, the 
traditional method, which simply excludes contaminat¬ 
ed pixels from the CMB map, does not fully eliminate 



the contribution of the point sources, which are mani¬ 
fest as small non-Gaussian features determined by the 
magnitude of the CMB fluctuations in the contaminat¬ 
ed locations. In addition, the sharp edges of the cutout 
sections lead to high-frequency noise that extends far 
beyond the frequency edge,which corresponds in our 
case to ^^aa,=1536. 

We illustrate the effect of applying our method in 
pure form in Fig.4, which shows (a) the angular power 
spectra of the convolved map CMB * BEAM without 
(solid curve) and with holes, whose size is determined 
by the size of the antenna beam (dashed curve), and 
(b) the same spectra for the input CMB map. We can 
see from the curves that the presence of holes in the 
maps leads to errors of about 13%i n the angular pow¬ 
er spectrum of the CMB for multipoles ^=200-300 and 
errors of about 30% for multipoles ^=450-550. The rel¬ 
ative errors are even higher for higher £ values. 

Note that, in spite of the fact that the radio sources 
considered are point sources (i.e.,they are essentially 
S functions), we have taken the contaminated loca¬ 
tions to include a region surrounding the source coordi¬ 
nates,determined by the characteristics of the inverse 
filter,which does not provide an ideal reconstruction 
of a ^ function. The presence of even a small amount 
of input noise requires regularization of the algorithm, 
which leads to a solution with finite resolution. To en¬ 
sure reliable exclusion of all contaminated pixels,we 
have assumed that the region distorted by each point 
source occupies an area equal to the area of the base 
of the antenna beam, which exceeds somewhat the to¬ 
tal area of the excluded CMB measurements. However, 
this is not a serious concern for our method,since we 
reconstruct the intrinsic CMB values in the resulting 
holes. 

Obviously, the maximum effect from reconstructing 
the CMB signal in contaminated locations is achieved 
in the absence of pixel noise. The quantitative effect is 
lowered as the pixel-noise level is increased, although 
the qualitative effect of a complete elimination of non- 
Gaussian features from the CMB maps is maintained 
even in the presence of fairly high noise levels. 

Figures 5a and 5b present the angular power spec¬ 
tra of the signals corresponding to experiments 2.1 and 
2.4, respectively.The results obtained in experiments 
2.2 and 2.3 occupy intermediate positions and are not 
shown. The solid curve shows the angular power spec¬ 
trum of the initial CMB map,the dotted curve the mea¬ 
sured CMB map, which satisfies (3), the dashed curve 
the angular spectrum of the reconstructed CMB map 
with the pixels distorted by point sources zeroed,and 
the dash-dotted curve the CMB map with reconstruct¬ 
ed values for the fluctuations at the distorted locations. 

We can see from Fig.5a that it is possible to recon¬ 
struct the CMB map nearly perfectly in the absence 
of pixel noise. The error in the reconstruction is deter¬ 


mined only by the error in the inversion filtration. The 
effect of applying our method in this case is close to 
that shown in Fig.4. 

Figure 5b shows that a high level of pixel noise 
leads to a substantial loss in the accuracy of the re¬ 
constructed CMB angular power spectrum,for both a 
simple exclusion of contaminated sections and the re¬ 
construction of the CMB signal in holes.However, the 
accuracy is nonetheless higher in the latter case — quite 
appreciably for multipoles ^=24-500. The dispersion of 
the residual-noise map for the specified interval of the 
input pixel noise is 25-27% lower due to the full re¬ 
moval of point sources from the CMB maps (Table 2). 

Although the simulations show that the effect of ap¬ 
plying our method decreases with growth in the pixel- 
noise level in a natural way, it nonetheless remains fair¬ 
ly high for a system with noise characteristics similar 
to those planned for the PLANCK mission [11]. 

The method for reconstructing the CMB signal pro¬ 
posed here can also be applied in connection with con¬ 
structing catalogs of point sources. For this, it is suf¬ 
ficient to subtract the reconstructed CMB map from 
the estimated CMB + PS map. 

6.3. Reconstruction of the CMB Anisotropy 
over a Wide Region of the Map 

From the point of view of enhancing the accura¬ 
cy of an estimated angular power spectrum based on 
CMB measurements over the entire celestial sphere, 
it is of interest to consider the reconstruction of the 
CMB signal in the zone of the Galaxy where powerful 
non-Gaussian noise is observed. 

In the traditional strategy, this region of the sky is 
simply not taken into account, which clearly leads to 
a loss of accuracy in the derived CMB angular power 
spectrum compared to that which would be obtained 
using a full set of undistorted data. Let us consider 
the reconstruction of the CMB signal in the zone of 
the Galaxy using the same method as in the previous 
experiment. Let the noise occupy the middle part of 
the map in the form of a band that is elongated in 
the horizontal direction. We cut out the contaminated 
region from the map and reconstruct the absent com¬ 
ponents of the CMB in this region using our method 
(experiment 3). 

Our simulations show that it is possible to obtain 
a nearly perfect reconstruction of the CMB, right up 
to the case when the width of the contaminated band 
comprises a third of the linear size of the map. Further 
increase in the width of the band leads to a growth in 
the errors of the reconstruction. Simulation results for 
band widths comprising about 30% of the linear size 
of the map are presented in Fig.6, which shows the 
(a) initial CMB map, (b) CMB map with the Galactic 
band cut out, and (c) reconstructed CMB map. The 
angular power spectra are presented in Fig.7, where 



the solid curve corresponds to the initial CMB map, 
the dashed curve to the map with the band cut out, 
and the dash-dotted curve to the reconstructed map. 

Analysis of the maps presented in Fig.7 shows that 
simple elimination of sections of the CMB occupying 
a substantial area leads to appreciable errors in the 
derived angular power spectrum, which can reach 50% 
(for example, for multipoles ^=200-300). Applying our 
method led to a virtually perfect reconstruction. 

Thus, this experiment provides hope that recon¬ 
structing the CMB in the zone of the Galaxy based on 
real measurements over the entire celestial sphere can 
substantially improve the accuracy with which the real 
angular power spectrum of the CMB can be derived. 

7. CONCLUSION 

We have investigated the fundamental possibility 
of improving the accuracy of estimates of the angu¬ 
lar power spectrum of the CMB by reconstructing the 
CMB anisotropy in regions contaminated by point ra¬ 
dio sources and other high-multipole noise. A series 
of examples have been used to demonstrate the effect 
of applying this approach, compared to the current¬ 
ly standard strategy of simply excluding contaminated 
sections of CMB maps. 

We have shown that, in the absence of instrumen¬ 
tal noise,the proposed method ensures the removal of 
point sources and other non-Gaussian multipole noise 
with very high accuracy,and is fairly stable to varia¬ 
tions in the noise level in the input data. 

The main advantage of the method is that it is 
able to fully remove non-Gaussian features from CMB 
maps, while the simple exclusion of contaminated lo¬ 
cations does not achieve this goal. The full removal of 
residual point sources from the maps makes it possi¬ 
ble to appreciably lower the dispersion of the residual 
noise, even in the presence of a relatively high level of 
input pixel noise. Our simulations show that the effect 
of applying our method decreases in a natural way as 
the input noise level grows, but nonetheless remains 
fairly high for noise characteristics corresponding to 
those planned for the PLANCK mission. 

We have also shown the fundamental possibility 
of enhancing the accuracy of estimates of the angu¬ 
lar power spectrum of the CMB by reconstructing 
the CMB anisotropy in the zone of the Galaxy where 
the strongest level of background contamination is ob¬ 
served. 
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Table 1. Dispersion a (in units of 10 of the maps in experiment 1 


CMB 

PS 

CMB + PS 

CMB with holes 

Residual 

map 

CM Brecon 

Residual 

map 

C M Brecon 
SNRfa 10) 

Residual 

map 

4.10 

4.75 

6.23 

3.97 

0.99 

4.10 

0 

4.11 

0.31 


Table 2. Maps and their dispersions for experiment group 2 


Experiment 

Map 

Dispersion a, 10 ^ 

Dispersion of map with holes a, 10 ^ 

2 

CMB 

4.10 

3.90 


PS * BEAM 

1.99 

0 


CMB* BEAM 

4.51 

4.29 


{CMB + PS) * BEAM 

4.93 

4.29 

2.1 

Pixelnoise 

0 

- 


C M Bmeas 

4.93 

- 


C M Brecon 

4.09 

3.89 


Residual map 

0.39 

1.28 

2.2 

Pixelnoise 

0.65 

— 


CM Bmeas 

4.97 

- 


C M Brecon 

4.05 

3.86 


Residual map 

0.74 

1.38 

2.3 

Pixelnoise 

2.61 

— 


C M Bfneas 

5.58 

— 


C M Brecon 

4.27 

4.10 


Residual map 

1.19 

1.59 

2.4 

Pixelnoise 

5.23 

— 


C M Bmeas 

7.19 

- 


C M Brecon 

4.14 

3.93 


Residual map 

1.47 

1.84 
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Fig.l. Maps for experiment 1. 
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Fig. 2 . Angular power spectra for the maps for experiment 1. See text for details. 
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Fig. 6 . Maps for experiment 3. 
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Fig.7. Angular power spectra for the maps for experiment 3. See text for details. 
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